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Abstract. We consider the transformation properties of integer sequences arising 
from the normal ordering of exponentiated boson ([a, a^] = 1) monomials of the form 
exp[A(a^) r a], r — 1,2,..., under the composition of their exponential generating 
functions (egf). They turn out to be of Sheffer-type. We demonstrate that two 
key properties of these sequences remain preserved under substitutional composition: 
a) the property of being the solution of the Stieltjes moment problem; and b) the 
representation of these sequences through infinite series (Dobihski-type relations). We 
present a number of examples of such composition satisfying properties a) and b). 
We obtain new Dobihski-type formulas and solve the associated moment problem for 
several hierarchically defined combinatorial families of sequences. 
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1. Introduction 



In a recent series of articles [1], [2], [3], [4], [5], [6] we investigated the properties of integer 
sequences appearing in the process of the normal ordering of powers of boson monomials 
[(a^) r a s ] n , with n, r, s -integers, where a and a) are the boson annihilation and creation 
operators respectively, satisfying [a, a)] = 1. They are extensions of earlier works [7], [8]. 
We observed that the normal form of \{a)) r a s ] n , with all the annihilation operators to 
the right, denoted by J\f ([(a t ) r a s ] n ) , can be written in the form (r > s): 

ns 

[(a r ) r a s ] n = M ([(a r ) r a s ] n ) = (a t ) n(r " s) J2 S rA n , k){a^) k a k i 1 ) 

k=s 

where S T:S (n, k) are generalizations of the conventional (r = s = 1) Stirling numbers of 
the second kind and 

ns 

B r , s (n) = Y,SrAn,k) (2) 

fc=s 

generalize the conventional (r = s = 1) Bell numbers. 

For general r > s we have worked out a complete theory of the numbers S rtS (n, k) 
and B ryS {n), including their recurrence relations, generating functions and closed- form 
formulas. In particular, the generalized Bell numbers B r;S (n) can be expressed as infinite 
series, thereby extending the celebrated Dobihski relation valid for r = s = 1 [9] : 

5i,i(») = -5]^ n = 0,l,2,... (3) 

fc=0 

Here are some examples of such relations: 

1 oo 1 n 

5 r> i(n) = -X; iT n[fc + (j-l)(r-l)] (4) 



fc=l j=l 



(5) 



fc=0 

they are all derived from the general polynomial-type formula (n — 1, 2, . . .) 



B r , s (n,y)= ^ S r . :S (n, k)y h 



k=s 

oo 



= e_?/ E^n[( A; +^- i )( r - s ))-( A; +^- i )( r - s )- 1 ) • 



fc! . 



...■(k + (j-l)(r-s)-s + l)}y k . (6) 
We may associate a Generating Function C(x) with a given sequence {c n } by [9] 



^) = E c «^- ( 7 ) 



?1=0 



This particular form of Generating Function is known as a Generating Function of 
Exponential Type or egf for short, due to the n\ denominators. Of particular interest for 
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us here are those sequences {B rs (n)} for which the egf can in fact be expressed as an 
exponential function; they include B r l (n), r = 1,2, ... for which 



e-^Y.BM^ (8) 

and [1],[2],[10] 



n=0 



cxp 



= -l) =Y B rl {n) — . r = 2,3,... (9) 

The numbers S r ^ s (n, k) appear when in Eqs.(8) and (9) an indeterminate y is introduced 
through 



^ 1l EE 5 mK%^ (io) 

n=0 \fc=l 

and 



cxp 



(\ "I oo / n \ n 

-E EwvK, (Hi 
n=0 \fc=l / 



2,3, 



Eqs.(10) and (11) define polynomials of order n: 

n 

BrM V) = Y1 S ^( n > k )y K r = 1, 2, . . . . (12) 
fc=l 

Evidently, B r ^(n) = B r . : i(n, 1). The polynomials of Eqs.(12) share another characteristic 
property: they can be written as ratios of two infinite series in y. These are the so-called 
Dobihski-type relations [1],[2], which for r = 1 and r > 1 respectively are: 

^Efc!^ = n = 0,l,... (13) 

k=l ' k=l 

and 

(t \\ n 00 T(n + —) n 

^E^frf^E^OaV, » = i,2,.... (i4) 

k=l ^^yr-l) k =l 

By setting y = 1 in Eqs.(10) and (11) we obtain a representation of the integers B r ^(n) 
as an infinite series (compare Eqs.(3)-(5)); this constitutes a fertile ground for their 
probabilistic interpretation [11], [12]. The numbers B r;1 (n) can also be given various 
combinatorial intepretations [13], [14]. The second consequence of Eqs.(13) and (14) 
(and of the more general formulas for s > 1, see [1],[2] ) is the fact, that B r ^(n,y) for 
y > is the n-th Stieltjes moment of a non-negative probability distribution, which is 
either discrete (for r = 1, giving a so called Dirac comb [3]) or continuous (for r > 1). 
This fact permits one to use the B rjS (n, y) to construct various quantum collective states 
called coherent states [4], [15]. The interpretation of combinatorial sequences as moments 
[5] has led to new calculational approaches to hyperdeterminants [16]. Another aspect of 
Eqs.(13) and (14) which deserves mention here is that the numbers SV,i(n, k) (1 < k < n) 
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form a non-singular lower-triangular matrix with ones on the diagonal. Such matrices 
form a group, called the Riordan group, which has important applications in enumerative 
combinatorics [17], [18]. 

The purpose of this note is to place Eqs.(10)-(14) in the more general context of 
Sheffer-type polynomials and to address the question of compositional substitution and 
its implication for the existence of Dobihski-type relations as solutions of the Stieltjes 
moment problem. 

We first recall the known fact [19], [20], [21] that a compositional substitution 
corresponds to multiplication of the matrices SV,i(n, k). Then we go on to demonstrate 
that if two polynomial sequences Bp(n,y) and Bc(n,y) generated by e yF ^ and e yG ^ 
respectively are solutions of the associated Stieltjes moment problems, then the sequence 
BF(G){n,y) is also a solution of another, closely related, Stieltjes moment problem. We 
further prove that if Bp(n, y) and Bc{n, y) are both given by Dobihski-type relations, see 
Eqs.(13) and (14), then the sequence B F ^(n,y) is also given by an analogous formula. 
We then illustrate these reproducing properties of Dobihski-type relations and moment 
problem solutions by some specific examples. They comprise multiple compositions of 
standard Bell numbers with themselves (composing discrete with discrete distributions), 
compositions of Lah numbers (related to Laguerre polynomials) with themselves and 
finally composing discrete with continuous distributions and vice versa. 

2. Sheffer-type polynomials 

A polynomial B(n, y) of order n in the variable y is of Sheffer-type if the associated egf 
can be written in the form [22] 



with A(0) = 1 and F(0) = 0. Many such families of polynomials have been thoroughly 
investigated. Among the polynomials encountered in Quantum Mechanics, the Hermite 
and Laguerre polynomials are of Sheffer-type, whereas the Legendre and Gegenbauer are 
not. Comparing Eq.(15) with Eqs.(10) and (11) we observe that B r> i(n,y) are Sheffer- 
type polynomials with A(x) — 1. In fact B 1: i(n, y) are the so-called Bell (or exponential) 
polynomials [22] and B 2 ^{n,y) are the generalized Laguerre polynomials. The numbers 
S ^(n, k) are the conventional Stirling numbers of the second kind and the numbers 



are the so-called unsigned Lah numbers [1],[10]. 

More generally, consider two families of Sheffer-type polynomials B F (n,y) and 
Ba(n,y) generated by 




(15) 



n=l 




(16) 




(17) 
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and 



oo / n 



= l + £ X>(n,% fc Up (18) 

n=l \fe=l / 

respectively, where F(0) = G(0) = and 

n n 

B F (n,y) = J2S F (n,k)y k and B G (n, y) = ^ 5 G (n, k)y k . (19) 
fc=i fc=i 

We now consider the polynomials generated by F(G(x)), i.e. 



oo In 



ey^ = l + J2[j:S F{G) (n,k)y k ) X - (20) 

n=l \fc=l / 

Before we calculate this sum we note the relation resulting from the change of summation 
in Eq.(18): 

oo / oo 



e yG(x) 



!+E E^-^ ^ ( 21 ) 

k=l \n=k J 

Now comparison with the direct expansion of the left hand side of Eq.(21) 

oo 

e yG(x) = i + J2(G(x)) k y k /k\ (22) 



k=i 

yields 

n=fc 

Proceeding to the direct calculation of Eq.(20) we recall that the matrices S F (n, k) and 
Sc(n, k) are lower triangular (i.e. the entries for k > n are zero): 



n=\ \k=l 
oo / n 



1 + E IE Sf ^ k ^ k E s ^ ^ (25) 

n=l \fc=l / p=n 



= 1 + E E E s c& n ) s ^ k ))y k ) -(/ ( 26 ) 

p=l \k=l \n=l J / P ' 

Comparison with Eq.(20) yields 

n 

S F ( G )(n, k) = ^ S G (n,p)S F (p, k) (27) 
P =i 

This last equality means that compositional substitution within the Sheffer-type 
polynomial families is equivalent to the matrix product of the corresponding Stirling 
matrices [19], [20] 

§F(G) = §>G • §F • (28) 
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A direct consequence of Eq. (27) is the formula 



n n 



B F(G) (n,y) = J2SF( G )(n,p)y p = J2y P J2 S ^ k ) S ^p) (29) 

p=l p=l k=l 

n k n 

= J2s G (n,k)J2SF(k,p)y p = J2s G (n,k)B F (k,y). (30) 

k=l p=l k=l 

The last equation can be seen as the generalized Stirling transform [23] of the 
polynomials B F (k,y) which for y — 1 reduces to the generalized Stirling transform 
of the sequence B F (k): 

n 

B F(G)(n) = S c(n, k)B F (k). (31) 
k=i 

3. Compositional moment problem 

The formulas (27) and (30) lead to important consequences if the initial Sheffer-type 
polynomials are solutions of the Stieltjes moment problems, i.e. if for x, y > there 
exist positive weight functions W F (x,y) and Wc(x,y) such that 

POO 

B F (n,y) = / x n W F (x,y)dx, (32) 
Jo 

POO 

B G (n,y)= / x n W G (x,y)dx. (33) 
Jo 

Then the following equalities follow: 

n 

B F{G) (n,y) = ^Scin, k)B F (k,y) 

k=l 

n POO POO n 

= ^2S G (n,k) / x k W F (x,y)dx = / W F (x, y) ^ S G (n, k)x k dx 
k=i Jo Jo k=i 

poo poo poo 

= / W F (x,y)B G (n,x)dx = / dxW F (x,y) / z n W G (z,x)dz 
Jo Jo Jo 

= So Zn (/ WF ( x > y ^ WG ( z > x ^ dx ) dz ( 34 ) 
and this implies that 

POO 

B F (G)(n,y)= x n W F(G) ( y x,y)dx (35) 
Jo 

where W F ( G )(x,y) is a positive function given by 

POO 

W F{G) (x,y)= W F (z,y)W G (x,z)dz. (36) 
Jo 

We remark that the arguments of the weight functions in Eq.(36) need not satisfy any 
particular symmetry properties. 



Hierarchical Dobinski-type relations via substitution and the moment problem 



7 



More generally, for p-fold substitution F 1 (F 2 (. ■ ■ (F p ) . . .)) one obtains 

/»oo /»00 

W Fl ( F2 (...(F p )...))(x,y) = / dz 1 W Fl (z 1 ,y) dz 2 Wp 2 (z 2) z 1 ) . . . 

Jo Jo 

poo 

... j dz p W Fv _ x {z P , z p _^)W Fv {x, z p ). (37) 
Jo 

Eq.(37) reveals a typical structure appearing in the iterated-kernel method of solving 
integral equations [24], [25]. 

In other words; for the Sheffer-type polynomials the property of being a solution 
of the Stieltjes moment problem is reproduced by the mechanism of compositional 
substitution, under the evident condition that the integrals in Eqs.(36) and (37) exist. 
In the following section we provide a number of examples of substitutions F(G(x)) for 
which an explicit evaluation of W F (G)( x ,y) and B F ^(n,y) can be carried through. 

4. Compositional Dobinski-type relations 

A rather large reservoir of solutions of the Stieltjes moment problem is contained in 
the formulas (13) and (14). For any r = 1,2, .. . B rt i(n,y) is the moment of a positive 
function W r (x,y), which can be written down explicitly, for instance by extending to 
y ^ I the results given in [4], [5], [6]. The examples are: 

W l{x ,y) = e-'±t^zR (38) 



k\ 

k=l 



W 2 {x,y) = ye 1 ^ y> j= — 

V^y 

W 3 (x, y) = ^—-e-%-v y (6V2x~ + SxyV^ F 2 & 2; 



(39) 



+v / 2x 3 /Vi^(l;|,2,|;^)), (40) 



where 8{z) is the Dirac delta function, I„(z) is the modified Bessel function of first 
kind and 0-^2 and 1F3 are hypergeometric functions. Eqs.(39) and (40) were obtained 
using the inverse Mellin transform. See [26] for its exposition and [27] for examples of 
applications. 

Note, that whereas Wi(x, y) is a discrete distribution in the form of a Dirac comb 
concentrated on positive integers, the functions W r (x,y) for r > 1 are continuous 
distributions [6]. Observe also that they are not normalized, in the sense of their zero 
moments: J °° W±(x, l)dx = 1 whereas J °° W r (x, l)dx 7^ 1, r > 1. 

In this section we demonstrate that the reproducing character of the compositional 
moment problem, see Eq.(36), implies the reproducing character of the Dobinski-type 
relations. In the following paragraph, with given F{x) and G(x) of Eqs.(17) and (18) we 
will carry out explicit substitutions F(G(x)) and analyze the weight functions W F (g)(x) 
obtained from Eq.(36) and the resulting Dobinski-type relations. 
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4.1. F(x) = G(x) = e x -I 

In the following the subscript B(B) stands for "substitute Bell into Bell" . We investigate 
the polynomials B B ( B )(n,y) resulting from 

e^ e ^ 1 - 1 )=f;^)(n )2/ )^ (41) 

n=0 

which correspond to the ordinary Stirling transform [23] of the Bell polynomials 

Bi,i(n,y) 

n 

B B (B)(n,y) = J2S(n,k)B 1A (k,y) (42) 
k=i 

where S(n, k) are the conventional Stirling numbers of the second kind. The polynomial 
B\,i(n,y) is the n-th moment of the Dirac comb [3], 

WB(x , y) = e -,£y^l, ( 43) 

k=i 

and the weight function resulting from the substitution F(F(x)) is through Eq.(36) 
equal to 



/•oo 

W B ( B )(x,y)= / W B (z,y)W B (x,z)dz 
Jo 



00 



e~ y 



k=i ' /V P =i y ' 



p=l F ' \k=l 



00 r / \ / P 

o(x — p) 



= e .^-i)££^ ^ 5(p , r )( ye -i r , (44) 

p=l \r=l / 

where the last equality results from the original Dobihski formula Eq.(13). This result 
shows that 

00 , n I v \ 
S B(fl) (n) = S B(B) (n, 1) = e( e_ " 1} £ ^ £ ^ ' ( 45 ) 

fc=l ' \r=l / 

with the initial terms B B(B) (n) = 1,1,3,12,60,358,2471,19302,..., for n = 0, 1, . . . . 
B B ( B ){n) counts the number of partitions of a set of n distinguishable elements, in which 
every part is again partitioned [19]. 

Multiple substitutions of Bell egf's into themselves result in hierarchical, chain- 
like formulas for corresponding partition numbers, i.e. for F(F(F(x))) one obtains for 
n = 0,l,... 



BB « B) ,(n)=e(^--' , -l)f;|. 



k=l 

k 



S(k, p)e~* S & 0e r(e_1 - 1) ) ) • (46) 

v p=l \r=l 
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For example, B B{B(B)) (n) = 1,1,4,22,154,1304,12915,146115,..., for n = 0,1..., 
which counts the number of "triple" partitions of an n-set. 

We conclude that the substitution F(F(x)) results in a formula for B B ^{n) which 
conserves the original Dobihski-type structure of B B (n) as in Eq.(3); and also gives 
a Dirac comb type of weight function with modified weights concentrated on positive 
integers. These results also hold good for higher order substitutions. 

4.2. F(x) = G(x) = 

This case corresponds to B 2t i(n,y) which from Eq.(14) is 

and can be also written as 

B 2A (n,y) = (n-l)\yL^U-y) (48) 

by using the standard form of the generating function of generalized Laguerre 
polynomials L^\x). With the notational convention introduced above we rewrite 
Eq.(48) as (here L stands for Laguerre) 

n 

B L (n,y) = J2S L (n,k)y k (49) 

k=l 

where Siin^k) are the unsigned Lah numbers, see Eq.(16). For y — 1, the integers 
B L (n, 1) = B L (n) count binary ordered forests of n nodes [13] (the initial terms are 
Bl(ji) = 1, 3, 13, 73, 501, 4051 . . . , n — 1, 2, . . .). For other combinatorial interpretations 
see [28]. 

The polynomial B L (n,y) is the n-th moment of [6] (see Eq.(40) ): 

W L (x,y)=ye-^ I -^0^- (50) 

v x y 

By F(F(x))-type composition the function exp (j^) generates B L ( L )(n,y) through 



<?^ = Y, B m(n>y)^ ( 51 ) 

n=0 

where L(L) stands for "substitute Laguerre into Laguerre" , which are the n-th moments 
of 



POD 

W L{L) (x,y)= / W L (z,y)W B (x,z)dz 
Jo 



oo 



<z+y) h(^/zy)\ ( _ (x+ , ) /i(2V5i) 



ye- (z+v > v ■ ze- {x+z > ±y 1Z ' dz. (52) 



By virtue of the entry 2.15.20.8 of [29], this yields a continuous distribution 

W L[L) (x,y) = ye~^^§± = \W L {% !), (53) 
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thus preserving the original structure encountered in Eq.(50). In addition, simple use 
of the generating function of the generalized Laguerre polynomials yields 

POO 

B L(L) {n,y) = / x n W L(L) {x ) y)dx = 2 n -\n-l)\yL^ l {-y i ) (54) 
Jo 

whose initial terms for y = 1 are -Bl(l) (n) — 1) 5, 37, 361, 4361, 62701 . . ., n — 1, 2, . . . . 
The p-fold substitution, p — 1,2, ... , gives in this case the compact expression: 

SL(L(...(L)...))(n) = V n ~\n - ly.L^U- 1 -), n = 1, 2, . . . (55) 
4.3. F(x)=e x -1, G(x) = j^ 

Here we substitute Laguerre (continuous distribution) into Bell (discrete distribution) 
and vice versa. 

The calculations are analogous to those in 4.1 and 4.2 with repeated use of integrals 
listed in [29]. We only quote the final results: 



B B 

where 



poo 

( L )(n,y)= x n W B {L){x,y)dx, (56) 
Jo 



W B(L) (x,y) = 1— J]| F> /fc e -*/ 1 (2v^) (57) 
^ x k=i 

which is a continuous distribution. The polynomials B B ^{n,y) are generated by 

oo „ 

e y^-i)^B B{L) {n,y)-. (58) 

n=0 

The initial terms of B B{L) {n) are 1, 4, 23, 171, 1552, 16583 . . . for n = 1, 2, . . .. These 
integers count structures called sets of sets of lists, where list means an ordered subset 
[28]. A closed-form Dobihski-type formula for B B ^(n) can be obtained by calculating 
the moments of W B ^(x, 1). A longer but straightforward calculation gives 

B B{L ){n) = e ^ (k^T)\ • (59) 

fe=i ^ >' 

Higher order substitutions yield formulas of similar type. 

For the opposite substitution ("Bell into Laguerre" denoted L(B) ) generated by 

e y 2=* = Y,B L(B) (n,y)? 1 . (60) 



TV. 

n=0 



we obtain 



where 



poo 

B L(B) (n,y)= / x n W L{B) (x,y)dx, (61) 
Jo 

KW*.»> = H£^r L "<-» (62) 

fc=i 
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which is a discrete (Dirac comb) distribution, with moments 

I 00 h.n-1 

B L{B) {n) = B L(B) (n, 1) = -e"* £ —Lg^-J) (63) 

fe=i 

and initial terms S B(L) (n) = 1, 4, 23, 173, 1602, 17575 . . ., for n = 1, 2, . . . . 

5e/Z numbers vs. "ordered" Bell numbers 

As the last example we shall consider a slightly more general substitution problem in 
which only the "internal" egf G(x) is of Sheffer-type. In other words, the egf of one of 
the sequences is not an exponential. A case in point is given by the so called "ordered" 
Bell numbers [3], [9] B (n) defined through 



B {n) = ^2S(n,k) k\ . (64) 
fc=i 

Their extension to polynomials B (n,y) = J2k=i ^( n ' ^) ^ V k * s generated by [21] 

1 00 n 

y v ' n=0 

Thus the B (n,y) are not of Sheffer-type. 

We now perform the substitution "Bell into ordered Bell" , denoted by the subscript 
0(B). Although Eq.(30) is no longer valid, we can still define the numbers B (B)(n) 
through Eq.(31): 

n 

B 0{ B)(n) = J2 S ( n > k ) B o(k), (66) 
k=i 

or equivalently by 

1 00 n 

^_ = J2B 0IB) (nf (67) 

ra=0 

Recalling the Dobihski-type expression for Bo(n) [3], [9] 

B o{n) = ~ 2 Y. Yk , (68) 

fc=0 

the formula Eq.(36), now for y = 1 only, carries over and after straightforward 
calculation we obtain the Dobihski-type formula for Bo(B)(n): 

fc=0 

where Li m (y) is the polylogarithm of order m of y. The initial terms are Bo(B){n) = 
1, 4, 23, 175, 1662, 18937, . . ., n = 1, 2, . . . . 

Similarly, from the substitution "double Bell into ordered Bell" (denoted by 
0(B(B)) below) we obtain 

50(B( B ))(n) = - £ If E Li -rd) ( 7Q ) 

fc=0 ' \r=l ' / 
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where 



1 



oc 



.71 



B o(B(B)){n) 



X 



(71) 



i 



I 



etc. 



2-e 



n=0 



Clearly, Eqs.(69) and (70) again give rise to Dirac comb weight functions. 
5. Discussion and conclusions 

The main result of this work can be viewed from different perspectives. It is primarily 
a method for the generation of new solutions of moment problems. As such it is of 
potential importance for the construction of new generalized coherent states. Refs. [5] 
and [4] should be considered as first steps in this direction. The iterative method based 
on Eqs.(36) and (37) appears to be straightforward under the condition of the existence 
of the relevant integrals. This will definitely extend and enrich the families of currently 
known solutions of the moment problem. 

A closer look at the examples above based on Eq.(36) leads to the conclusion that 
if e G( - x ^ generates the moments of a discrete distribution then the moments generated 
by e F ( G ( x ^ are those of a discrete distribution. Similarly, when e G ^ gives a continuous 
distribution, the composition e F( - G( - x ^ gives rise to a continuous distribution. 

We are dealing here with Sheffer-type polynomials which are also solutions of the 
moment problem; it should be borne in mind that these are quite strong restrictions. 
It is easy to construct Sheffer-type polynomials which are not solutions of the moment 
problem. For example, the polynomials p n (y), which are related to Bessel polynomials 
[22], are generated by 



and can take on negative values for y — 1; they are therefore not acceptable solutions of 
the moment problem. On the other hand, for s > 1, the polynomials B r s (n,y) defined 
by Eq.(6) are solutions of the moment problem [6] but are not of Sheffer-type [1],[2]. 

Referring to various Dirac comb-type distributions obtained by compositions (see 
Eqs.(43), (45), (46), (63), (68), (69) and (70) ) we observe that the substitution B(n) -> 
B(an 2 + j3n + 7), (a, /3, 7 - integers, a > 0) gives sequences B(n) = B(an 2 + f3n + 7) 
which are the n-th moments of continuous measures; they are infinite, weighted sums 
of log-normal distributions [3] . 

The reproducing nature of Dobihski-type relations under composition also follows 
from the scheme presented here. It has already provided a number of new closed-form 
expressions for combinatorial numbers, Eqs. (54), (55), (59), (63), (69) and (70), together 
with the associated weight functions. It seems that this method can be also applied to 
various generalizations of combinatorial numbers, e.g. q- deformations [30] and to more 
involved substitution schemes such as those considered in [31]. 




(72) 



n=l 
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